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Abstract 

Our main theoretical result is that, if a simple polytope has a pair of complementary 
vertices (i.e., two vertices with no facets in common), then it has a second such pair. 
Using this result, we improve adjacency testing for vertices in both simple and non-simple 
polytopes: given a polytope in the standard form {x £ R" | Ax — b and x > 0} and a list 
of its V vertices, we describe an 0{n) test to identify whether any two given vertices 
are adjacent. For simple polytopes this test is perfect; for non-simple polytopes it may 
be indeterminate, and instead acts as a filter to identify non-adjacent pairs. Our test 
requires an 0{n^V + nV^) precomputation, which is acceptable in settings such as all-pairs 
adjacency testing. These results improve upon the more general 0{nV) combinatorial and 
O(n^) algebraic adjacency tests from the literature. 
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1 Introduction 

Two vertices of a polytope are complementary if they do not belong to a common facet. Com- 
plementary vertices play an important role in the theory of polytopes; for instance, they provide 
the setting for the d-step conjecture [Hllin] (now recently disproved [H]), and in the dual setting 
of disjoint facets they play a role in the classification of compact hyperbolic Coxeter polytopes 
[B] . In game theory, Nash equilibria of bimatrix games are described by an analogous concept 
of complementary vertices in pairs of polytopes [5J [TB] . 

Our first main contribution, presented in Section [2l relates to the minimal number of com- 
plementary vertex pairs. Many polytopes have no pairs of complementary vertices at all (for 
instance, any neighbourly polytope). However, we prove here that if a simple polytope P of 
dimension d > 1 has at least one pair of complementary vertices, then it must have at least two 
such pairs. 

The proof involves the construction of paths through a graph whose nodes represent pairs of 
complementary or "almost complementary" vertices of P. In this sense it is reminiscent of the 
Lemke-Howson algorithm for constructing Nash equilibria in bimatrix games [T^] , although our 
proof operates in a less well-controlled setting. We discuss this relationship further in Section|31 

Our second main contribution, presented in Section [31 is algorithmic: we use our first 
theorem to build a fast adjacency test. Specifically, given a polytope in the standard form 



P = {x £ M" I = b and x > 0} with V vertices, we begin with an 0{n^V + nV"^) time pre- 
computation step, after which we can test any two vertices for adjacency in 0{n) time. If P 
is simple (which the algorithm can also identify) then this test always gives a precise response; 
otherwise it may be indeterminate but it can still assist in identifying non-adjacent pairs. 

The key idea is, for each pair of vertices u,v & P, to compute the join w V u; that is, the 
minimal face containing both u and v. If P is simple then our theorem on complementary 
vertices shows that m V u — u' V v' for a second pair of vertices u',v'. Our algorithm then 
identifies such "duplicate" joins. 

Although the precomputation is significant, if we are testing all (^) pairs of vertices for 
adjacency then it does not increase the overall time complexity. Our 0{n) test then becomes 
extremely fast, outperforming the standard 0{nV) combinatorial and O(n^) algebraic tests from 
the literature [7]. Even in the non-simple setting, our test can be used as a fast pre-filter to 
identify non-adjacent pairs of vertices, before running the more expensive standard tests on 
those pairs that remain. 

In Section |4] we discuss these performance issues further, as well as the application of these 
ideas to the key problem of polytope vertex enumeration. 

All time complexities are measured using the arithmetic model of computation, where we 
treat each arithmetical operation as constant-time. 

We briefly remind the reader of the necessary terminology. Following Ziegler |17| , we insist 
that all polytopes be bounded. A facet of a d-dimensional polytope P is a (d — l)-dimensional 
face, and two vertices of P are adjacent if they are joined by an edge. P is simple if every vertex 
belongs to precisely d facets (i.e., every vertex figure is a (d — l)-simplex), and P is simplicial 
if every facet contains precisely d vertices (i.e., every facet is a (d — l)-simplex). As before, two 
vertices of P are complementary if they do not belong to a common facet; similarly, two facets 
of P are disjoint if they do not contain a common vertex. 

2 Complementary Vertices 

In this section we prove the main theoretical result of this paper: 

Theorem 1. Let P be a simple polytope of dimension d > 1. If P has a pair of complementary 
vertices, then P has at least two distinct pairs of complementary vertices. 

Note that these pairs of vertices are distinct, but need not be disjoint: they might be of the 
form {u,v} and {u,w} for some vertices u,v,w € P. 

The proof of Theorem [T] involves an auxiliary graph, which we now describe. To avoid 
confusion with vertices and edges of polytopes, we describe graphs in terms of nodes and arcs. 
We do not allow graphs to have loops or multiple edges. 

Definition 2 (Auxiliary graph). Let P be a simple polytope of dimension d > 1. We construct 
the auxiliary graph r(P) as follows: 

• The nodes of r(P) are unordered pairs of vertices {u, v} of P where u, v have at most one 
facet in common. We say the node {u,v} is of type A if the vertices u,v G P are comple- 
mentary (they have no facets in common) , or of type B otherwise (they have precisely one 
facet in common). 

• The arcs of r(P) join nodes of the form {u,x} and {u,?/}, where x and y are adjacent 
vertices of P, and where no single facet of P contains all three vertices u, x, y. 
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Figure 1: A polytope P and the corresponding auxiliary graph T{P) 



Figure [T] illustrates this graph for the case where P is a cube. Informally, each arc of r(P) 
modifies a node by "moving" one of its two vertices along an edge of P, so that if the two 
vertices lie on a common facet F then this movement is away from F. More formally, we can 
characterise the arcs of r(P) as follows: 

Lemma 3. Let v = he a node oJT{P) as outlined above. 

(i) If v is of type A, there are precisely 2d arcs meeting v. These include d arcs that connect 
V with {u,x} for every vertex x adjacent to v in P, and d arcs that connect v with {y,v} 
for every vertex y adjacent to u in P. 

(ii) If V is of type B, there are precisely two arcs meeting v. Let F he the unique facet of P 
containing hoth u and v. Then these two arcs join v with {u,x} and {y,v}, where x is the 
unique vertex of P adjacent to v for which x ^ F, and y is the unique vertex of P adjacent 
to u for which y ^ F. 

Proof. We consider the type A and B cases in turn. 

(i) Let = {u,v} be of type A, and let x be any vertex of P adjacent to v. Since v is of 
type A, u and v have no facets in common. Since P is simple, at most one facet of P 
contains x but not v. Therefore u and x have at most one facet in common, and so {w, x} 
is a node of r(P). 

Since {u, x} is a node of r(P) and no facet contains both u and v, it follows that — {u, v} 
and {u, x} are joined by an arc. A similar argument applies to nodes {y, v} where y is any 
vertex of P adjacent to u. Because P is simple there are precisely d vertices adjacent to v 
and d vertices adjacent to u, yielding precisely 2d arcs of this type. 

(ii) Now let = {u, v} be of type B, and let F be the unique facet of P containing both u and 
V. If there is an arc from v to any node of the form {u,x}, it is clear from Definition [2] 
that we must have x adjacent to v and x ^ F. Because P is simple, there is precisely one 
X with these properties. 

We now show that an arc from v to {u, x} does indeed exist. Because P is simple, at most 
one facet of P contains x but not v. The vertex v in turn has only the facet F in common 
with It; since x ^ F it follows that u and x have at most one facet in common. Therefore 
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{u, x} is a node of T{P). Because x ^ F there is no facet containing all of u, v and x, and 
so the arc from v to {u,x} exists. 

A similar argument applies to the arc from ly to {y,v}, yielding precisely two arcs that 
meet v as described in the lemma statement. 

To finish, we note that every vertex belongs to d > 1 facets, and so every node {u, v} of 
r(P) has u ^ V. This ensures that we do not double-count arcs in our argument; that is, the 
2d arcs in case (i) join i' to 2d distinct nodes of r(P), and likewise the two arcs in case (ii) join 
1/ to two distinct nodes of T{P). □ 

Our strategy for proving Theorem [1] is to show that, if we follow any path from a type A 
node of r(P), we must arrive at some different type A node; that is, we obtain a new pair of 
complementary vertices. The details are as follows. 

Proof of Theorem]^ To establish Theorem [TJ we must prove that if T{P) contains at least one 
type A node, then it contains at least two type A nodes. 

Let v he a type A node, and let a be any arc meeting v. Since every type B node has 
degree two (by Lemma [3]), this arc a begins a well-defined path through r(P) that passes 
through zero or more type B nodes in sequence, until either (a) it arrives at a new type A node 
z/', or (b) it returns to the original type A node v and becomes a cycle (see Figure [2|). Case (a) 
gives us the desired result; our task is to prove that case (b) is impossible. 




O Type A 
O Type B 



Case (a) Case (b) 



Figure 2: Following a path from the type A node ly 

Suppose then that case (b) occurs, and there is a cycle that passes through nodes i^, /^i, /X2, . . . , 
Hk, V in turn, where v is type A and ^i, . . . , ^j,k are type B. Since r(P) is a graph with no loops 
or multiple edges, we have k > 2. Since deg(/ii) = 2, the nodes /ii, . . . , /i^ are all distinct. 

For any node rj = {u, v} of T{G), let $(77) denote the set of all facets of P that contain either 
u or V. Because P is simple, Definition [2] gives |<i>(?7)| = 2d if 77 is type A, or |<i>(77)| = 2(i — 1 if 
77 is type B. 

We now show that $(/ii) = ^■(^2) = ■ ■ ■ — ^{nk)- Consider two adjacent nodes /i^ = {x,y} 
and Hi+i = {x, z} along our cycle. Because y and z are adjacent vertices of the simple polytope 
P, there is only one facet F for which y E F but z ^ F. By Lemma [3l this facet F must be the 
unique facet containing both x and y. Therefore every facet that contains either x or y must 
also contain either x or z, and we have ^(^i) C $(/ij_|_i). Because [^(/ii)! — |$(/Xi+i)| = 2d — 1 
we have ^{fJ-i) = $(/Xi+i), and by induction $(^^1) — *&(M2) — . . . — ^{fJ,k)- 

Consider now the type A node ly that begins and ends the cycle. Let ly = {u, v}, and without 
loss of generality let = {u, x} for some vertex a; e P. Because |$(/xi)| = 2d—l < 2d — |$(i^)|, 
there must be some facet F for which F G $(7^) but F ^ ^(/ii); it follows then that v £ F but 
x^F. 
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Moving to the other end of the cycle: since is adjacent to v in r(P), we have either 
Mfc — {u, 2;'} where x' is adjacent to v in P, or else /ifc = {x' , where x' is adjacent to u in P. 
The second option is not possible because v € F ^ 'i'(Mi) = ^il^k)] therefore — {w, a;'} for 
some x' adjacent to v. 

We now know that both vertices x and x' are adjacent to v\ moreover, since F ^ ^{^l^i) = 
we have x, x' ^ F . However, P is a simple polytope with v G F, and so there is only one 
vertex adjacent to v that is not in F. Therefore x = x' and /ii = /i^, contradicting the earlier 
observations that k > 2 and the nodes /ii, . . . , /i^ are all distinct. This completes the proof of 
Theorem [H □ 

Remark. The proof of Theorem [T] is algorithmic: given a simple polytope P of dimension 
d > 1 and a pair of complementary vertices u,v & P, it gives an explicit algorithm for locating 
a second pair of complementary vertices. 

In essence, we arbitrarily replace one of the vertices u with an adjacent vertex u', and then 
repeatedly adjust this pair of vertices according to Lemma [3] part (ii) until we once again reach 
a pair of complementary vertices. For each adjustment. Lemma [3] part (ii) gives two options 
(corresponding to the two arcs that meet a type B node); we always choose the option that 
leads us "forwards" to a new pair, and not "backwards" to the pair we had immediately before. 

The proof above ensures that we will eventually reach a complementary pair of vertices again, 
and that these will not be the same as the original pair u, v. 

Passing to the dual polytope, Theorem [T] gives us an immediate corollary: 

Corollary 4. Let P he a simplicial polytope of dimension d > 1. If P has a pair of disjoint 
facets, then P has at least two distinct pairs of disjoint facets. 

We finishing by showing that the "simple" and "simplicial" conditions are necessary in 
Theorem [1] and Corollary |4l 

Observation 5. The triangular bipyramid (Figure\^ left) is a non-simple polytope of dimension 
d = 3 with precisely one pair of complementary vertices ( the apexes at the top and bottom, shaded 
in the diagram). 

Its dual is the triangular prism (Figure\3[ right), which is a non- simplicial polytope of di- 
mension d = 3 with precisely one pair of disjoint facets (the triangles at the top and bottom, 
shaded in the diagram). 



These constructions are easily generalised. For instance, we can build a non-simple bipyramid 
over a neighbourly polytope: the two apexes form the unique pair of complementary vertices, 
and all other pairs of vertices are adjacent. Felikson and Tumarkin provide further examples in 
the dual setting [6j, involving non-simplicial Coxeter polytopes with precisely one pair of disjoint 
facets. 






Figure 3: A triangular bipyramid (left) and a triangular prism (right) 
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3 Adjacency Testing 



In this section we prove our main algorithmic result, which uses Theorem [T] to identify simple 
polytopes and test for adjacent vertices. Throughout this section we work with polytopes of the 
form 

P = {x e M" I ylx b and X > 0}, (1) 

where A is some n-column matrix. This form is standard in mathematical programming, and 
appears in key applications of vertex enumeration jJl [S]. Note that the dimension of P is 
not immediately clear from ([T|), although n — rank A gives an upper bound; likewise, it is not 
immediately clear whether P is simple. 

We recall some standard terminology from polytope theory: if F and G are faces of the 
polytope P, then the join P V G is the unique smallest-dimensional face that contains both P 
and G as subfaces. For example, recall the cube from Figure [U and consider the edges AD and 
DC. Their meet AD A DC is the vertex D (where they intersect), and their join AD V DC is 
the square facet ABCD (the smallest face containing both edges). 

Note that the join may be the entire polytope P (for instance, this happens if the faces P 
and G are complementary vertices). 

Our main algorithmic result is the following: 

Theorem 6. Consider any polytope P = {x G R" | Ax — h and x > 0}, and suppose we have a 
list of all vertices of P, with V vertices in total. Then, after a precomputation step that requires 
0{n'^V + nV^) time and 0{nV^) space: 

• we know immediately the dimension of P and whether or not P is .simple; 

• if P is simple then, given any two vertices u,v ^ P , we can test whether or not u and v 
are adjacent in 0{n) time; 

• if P is not simple then, given any two vertices u,v ^ P , we may still he able to identify 
that u and v are non- adjacent in 0{n) time. 

We discuss the importance and implications of this result in Section |4l in the meantime, we 
devote the remainder of this section to proving Theorem [6l 

Our overall strategy is to use our main theoretical result (Theorem [T|) to characterise and 
identify adjacent vertices. As a first step, we describe complementary vertices in terms of joins: 

Lemma 7. Let u and v be distinct vertices of a polytope P. Then u and v are complementary 
vertices if and only if uW v = P. 

Proof. If u V V ^ P then there is some facet P for which uV v C P; therefore u,v G P, and u 
and V cannot be complementary. 

On the other hand, if it V w = P then there is no facet P for which u,v G F (otherwise we 
would have u V w C P) . Therefore u and v are complementary. □ 

Using this lemma, we can now make a direct link between Theorem [1] and adjacency testing 
in polytopes: 

Theorem 8. Let P be a simple polytope, and let F be a face of P. Then F is an edge if and 
only if there is precisely one pair of distinct vertices u,v € P for which F = u\/ v. 

If P is any polytope (not necessarily simple), then the forward direction still holds: if F is 
an edge then there must be precisely one pair of distinct vertices u,v £ P for which F = uV v. 
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Proof. The forward direction is straightforward: let F be an edge of any (simple or non-simple) 
polytope P, and let the endpoints of F be the vertices u and v. It is clear that F = u\/ v, and 
because F contains no other vertices it cannot be expressed as the join x V y for any other pair 
of distinct vertices x, y. 

We now prove the reverse direction in the case where P is a simple polytope. Let F be a 
face of P, and suppose that F is not an edge. If dimP < 1 then P contains at most one vertex, 
and so there can be no pairs of distinct vertices m, v for which F — u\/ v. 

Otherwise dimP > 1; moreover, since P is a simple polytope then the face P is likewise 
simple when considered as a polytope of its own. Hence we can invoke Theorem [1] to finish the 
proof. If F ~ uV V for distinct vertices u, v, then Lemma [7] shows that u, v are complementary 
in the "sub-polytope" P. By Theorem [T] there must be another pair of complementary vertices 
{u' , v'} ^ {u, v} in P, and by Lemma[7] we have F = u' M v' as well. □ 

To make the join operation accessible to algorithms, we describe faces of polytopes using 
zero sets. Fukuda and Prodon [7j define zero sets for points in P; here we extend this concept 
to arbitrary faces. 

Definition 9 (Zero set). Consider any polytope of the form P = {x £ M" | Ax = b and x > 0}, 
and let P be any face of P. Then the zero set of P, denoted Z{F), is the set of coordinate 
positions that take the value zero throughout P. That is, Z{F) = {i \ Xi = for all x G P} C 
{1,2,. 

Zero sets are efficient for computation: each can be stored and manipulated using a bit- 
mask of size n, which for moderate problems {n < 64) involves just a single machine-native 
integer and fast bitwise CPU instructions. Joins, equality and subface testing all have natural 
representations using zero sets: 

Lemma 10. Let F and G be non-empty faces of the polytope P = {x e R" | Ax. = b and x > 0}. 
Then: 

• F QG if and only if Z{F) D Z{G); 

• F ^G if and only if Z{F) ^ Z{G); 

• the jom has zero set Z{F V G) = Z{F) n Z{G). 

Proof. For each i = 1, . . . , n, let iJ^ be the hyperplane Hi = {x e K" \xi — Q}. It is a standard 
result that every non-empty face of P is the intersection of P with some set of hyperplanes Hi , 
and conversely that any intersection of P with some set of hyperplanes Hi yields a (possibly 
empty) face of P. 

We claim that every non-empty face P C P satisfies F = P (1 ^Pl^g^^j^-j Hi^ ; that is, P is 
obtained by intersecting P with every Hi for which i G Z(F). 

To prove this claim, consider the face F' = P H (C\i^z(F) -^i) ■ definition of Z{F) it is 
clear that F C F' . li F ^ F' then it follows that P is a strict subface of P', and so there must 
be some additional hyperplane Hj with j ^ Z{F) for which F C F' O Hj . This in turn would 
imply that j G Z{F), a contradiction. Therefore F = F' . 

We now prove the individual statements of the lemma. If P C G then it is clear by definition 
of the zero set that Z{F) D Z{G). Conversely, if Z{F) 3 Z{G) then it follows from the claim 
above that P is the intersection of G with zero or more additional hyperplanes Hi, and so 
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F C G. This proves the first statement, and the second statement of the iemma now foliows 
immediately. 

We finish with the third statement. Since F, G C F V G, it follows from the first statement 
that Z{F V G) C Z{F) n Z{G). Suppose now that Z{F V G) ^ Z{F) D Z{G); that is, there is 
some i G Z{F) n Z{G) for which i ^ Z{F V G). Denote F' = n (F V G). Since i ^ Z{F V G) 
we see that F' is a strict subface of F V G; moreover, since i G Z{F) n Z{G) we have F,G C Hi 
and so F,G Q F'. Therefore F V G is not the smallest-dimensional face containing both F and 
G, contradicting the definition of join. □ 

We now define the main data structure that we build during our precomputation step in 
Theorem [6l 

Definition 11 (Join map). Consider any polytope of the form P = {x g E" | Ax = b and x > 0} 

The join map of P, denoted Jp, is a map of the form Jp : 2^^' - '"^ Z>o; that is, J^p maps 
subsets of {!,...,«} to non- negative integers. For each subset S C {1, . . . , n}, we define the 
image Jp{S) to be the number of pairs of distinct vertices {u, f } G P for which Z{u y v) — S. 

The following result reformulates Theorem [8] in terms of the join map, and follows immedi- 
ately from Theorem [8] and Lemma 1101 

Corollary 12. Let P he a simple polytope, and let u,v be vertices of P. Then u and v are 
adjacent if and only if Jp{Z(u) n Z{v)) = 1. 

If P is any polytope (not necessarily simple), then the forward direction still holds: if u and 
V are adjacent then we must have Jp(Z(u) fl ^(f)) — 1. 

As a further corollary, the join map can be used to identify precisely whether or not a 
polytope is simple: 

Corollary 13. Let P he any polytope of dimension d. Then P is simple if and only if, for every 
vertex u G P, there are precisely d other vertices u' £ P for which Jp{Z{u) H Z(u')) = 1. 

Proof. If P is simple then, for each vertex u Cz P, there are precisely d other vertices m' G P 
adjacent to u. By Corollary [T^ it follows that there are precisely d other vertices u' E P for 
which Jp(Z{u) n Z{u')) = 1. 

If P is non-simple then there is some vertex u E P that belongs to > d edges, and so there 
are > d other vertices u' G P adjacent to u. By Corollary [121 we have Jp{Z{u) r\Z{u')) = 1 for 
each of these adjacent vertices. □ 

We now show that the join map enjoys many of the properties required for the complexity 
bounds in Theorem [6] 

Lemma 14. Consider any polytope P = {x G M" | Ax = b and x > 0}, and suppose we have 
a list of all vertices of P, with V vertices in total. Then we can construct the join map Jp in 
OiriV'^) time and 0{nV^) space. Once it has been constructed, we can compute Jp{S) for any 
set S C {1, . . . , n} in 0{n) time. 

Proof. We store the join map J^p using a trie (also known as a prefix tree). This is a binary 
tree of height n. Each leaf node (at depth n) represents some set S" C {1, . . . , n}, and stores the 
corresponding image Sp{S). Each intermediate node at depth k < n supports two children: a 
"left child" beneath which every set S has k + 1 ^ S, and a "right child" beneath which every 
set S has A; + 1 G S". 
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We optimise our trie by only storing leaf nodes for sets S with Jp{S) > 1, and only storing 
intermediate nodes that have such leaves beneath them. Figure 0] illustrates the complete trie 
for an example map with n = 3. 




Figure 4: An example of a trie representing a join map for n = 3 

Tries are fast to use: for any set S C {1, . . . , n}, the operations of looking up the value of 
Jp{S), inserting a new value of Jp{S), or updating an existing value of Jp{S) are all Oin), 
since each operation requires us to follow a single path from the root down to level n (possibly 
inserting new nodes as we go). For further information on tries in general, see a standard text 
such as fm . 

It is clear now that we can construct Jp in 0{nV^) time: for each of the (^) pairs of vertices 
{u,w}, we construct the set S = Z{u\/v) = Z{u)r\Z{v) in 0{n) time (just test which coordinate 
positions are zero in both u and w), and then perform an 0{n) lookup for S in the trie. If S 
is present then we increment Jp{S)\ otherwise we perform an 0{n) insertion to store the new 
value Jp{S) = 1. 

It is also clear that the trie consumes at most 0{nV'^) space: because the construction 
involves 0{V'^) insertions we have 0{V'^) leaf nodes, and since the trie has depth n this gives 
0{nV^) nodes in total. 

Finally, for any set S C {1, . . . , n}, computing Jp{S) involves a simple lookup operation in 
the trie, which again requires 0{rL) time. □ 

We are now ready to complete the proof of our main result. Theorem [6l 

Proof of Theorem\^ The precomputation involves three stages: 

(i) building the join map Jp; 

(ii) computing dim P using a rank computation; 

(iii) iterating through the join map to determine whether P is simple. 
By Lemma [HI stage (i) requires 0{nV^) time and 0{nV^) space. 

To compute the dimension in stage (ii) we select an arbitrary vertex vq G P, and build a 
{V — 1) X n matrix AI whose rows are of the form v — vq for each vertex v ^ vq. It follows 
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that dimP — rank A/, and using standard Gaussian elimination we can compute this rank in 
0{r?V) time and O(nF) space. 

For stage (iii) we once again scan through all (^) pairs of vertices m, v, compute Z(u) n Z(v) 
for each in 0(n) time, and use an 0(n) lookup in J7p to test whether Jp[Z{u) fl Z{v)) — 1. 
We can thereby identify whether or not, for each vertex u ^ P, there are precisely dimP other 
vertices u' € P for which Jp{Z{u) n Z{u')) = 1; by Corollary [T3] this tells us whether or not P 
is simple. The total running time for this stage is 0{nV^). 

Summing the three stages together, we find that our precomputation step requires a total of 
0{n'^V + nV"^) time and 0{nV'^) space. 

Once this precomputation is complete, it is clear from stages (ii) and (iii) that we know 
immediately the dimension of P and whether or not P is simple. 

Consider now any two vertices u,v & P. As before, we can compute the set Z{u) n Z{v) 
in 0{n) time, and using Lemma [T3] we can evaluate Jp{Z{u) fl Z{v)) in 0{n) time. Now 
Corollary [T^l tells us what we need to know: if P is simple then u and v are adjacent if and only 
if Jp{Z{u) n Z{v)) = 1, and if P is non-simple but Jp{Z{u) n Z{v)) ^ 1 then we still identify 
that u and v are non-adjacent. □ 



4 Discussion 

As noted in the introduction, the proof of Theorem [T] is reminiscent of the Lemke-Howson 
algorithm for constructing Nash equilibria [T2]. The Lemke-Howson algorithm operates on a 
pair of simple polytopes P and Q (best response polytopes for a bimatrix game), each with 
precisely / = dimP -I- dimQ facets labelled 1, . . . ,/, and locates vertices u e P and u G Q 
whose incident facet labels combine to give the full set {1, . . . , /}. See [21 HI] for details. 

The Lemke-Howson algorithm can also be framed in terms of paths through a graph F, where 
it can be shown that these paths yield a 1-factorisation of F. One then obtains the corollary 
that the number of fully-labelled vertex pairs is even; in particular, because there is always a 
"trivial" pair (0, 0), there must be a second pair (which gives to a Nash equilibrium). 

In this paper our setting is less well controlled. We work with a single polytope P, which 
means we must avoid transforming the pair of vertices into the identical pair {v,u}. 

Moreover, P may have arbitrarily many facets, which makes the arcs of F(P) more difficult 
to categorise. In particular, we do not obtain any such 1-factorisation or parity results. To 
illustrate this. Figure [S] shows that both parities can occur: the cube (on the left) has four 
complementary pairs of vertices, whereas the cube with one truncated vertex (on the right) has 
nine. 




Figure 5: Simple polytopes with (i) four, and (ii) nine complementary vertex pairs 

Moving to Theorem [6l our 0{n) time adjacency test is the fastest we can hope for, since 
vertices require VL{n) space to store (a consequence of the fact that there may be exponentially 
many vertices [13]). In contrast, standard approaches to adjacency testing use either an 0{nV) 
"combinatorial test" (where we search for a third vertex w G u V ?;), or an 0{'n?) "algebraic 
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test" (where we use a rank computation to determine dim(M V v)). See [7] for details of these 
standard tests. 

If we are testing adjacency for ah pairs of vertices, our total running time including precom- 
putation comes to 0{n^V + nV^), as opposed to 0{nV^) or 0{n^V^) for the combinatorial and 
algebraic tests respectively. This is a significant improvement, given that V may be exponential 
in n. 

The main drawback of our test is that it only guarantees conclusive results for simple poly- 
topes. Nevertheless, it remains useful in the general case: it can detect when a polytope is 
simple, even if this is not clear from the initial representation {x G K." | Ax = b and x > 0}, and 
even in the non-simple setting it gives a fast filter for eliminating non-adjacent pairs of vertices. 

A key application of all-pairs adjacency testing is in the vertex enumeration problem: given 
a polytope in the form {x S M" | Ax = b and x > 0}, identify all of its vertices. This is a 
difficult problem, and it is still unknown whether there exists an algorithm polynomial in the 
combined input and output size. 

The two best-known algorithms are reverse search [TJ [2] and the double description method 
[71 [13]: each with their own advantages and drawbacks. The double description method, which 
features in several application areas such as multiobjective optimisation i4l and low-dimensional 
topology [5], inductively constructs a sequence of polytopes by adding constraints one at a time. 
Its major bottleneck is in identifying pairs of adjacent vertices in each intermediate polytope, 
and in this setting our fast all-pairs adjacency test can be of significant practical use. 
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